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We investigate, both experimentally and theoretically, possible routes towards Anderson-like lo- 
calization of Bose-Einstein condensates in disordered potentials. The dependence of this quantum 
interference effect on the nonlinear interactions and the shape of the disorder potential is investi- 
gated. Experiments with an optical lattice and a superimposed disordered potential reveal the lack 
of Anderson localization. A theoretical analysis shows that this absence is due to the large length 
scale of the disorder potential as well as its screening by the nonlinear interactions. Further analysis 
shows that incommensurable superlattices should allow for the observation of the cross-over from the 
nonlinear screening regime to the Anderson localized case within realistic experimental parameters. 



Disordered systems have played a central role in con- 
densed matter physics in the last 50 years. Recently, it 
was proposed that ultracold atomic gases may serve as 
a laboratory for disordered quantum systems 0, Q and 
allow for the experimental investigation of various open 
problems in that field 0. Some of these problems con- 
cern strongly correlated systems the realization of 
Bose 0, or Fermi glasses , quantum spin glasses 13 
and quantum percolation |^ . This letter addresses one of 
the most important issues, namely the interplay of An- 
derson localization (AL) [lOj and repulsive interactions 
This interplay may lead to the creation of delocal- 
ized phases both for fermions and bosons 0. The 
possible occurrence of AL has also been investigated the- 
oretically for weakly interacting Bose-Einstein conden- 
sates (BEC) ^3 J an.d in this case it was shown that even 
moderate nonlinear interaction counteracts the localiza- 
tion. As a main result of this letter we show that despite 
this difficulty there exists an experimentally accessible 
regime where Anderson-like localization can be realized 
with present day techniques. 

Several methods have been proposed to produce a dis- 
ordered, or quasi-disordered potential for trapped atomic 
gases. They include the use of speckle radiation 0|, in- 
commensurable optical lattices [isj . impurity atoms in 
the sample and the disorder that appears close to 
the surface of atom chips Recently, first experi- 

ments searching for effects of disorder in the dynamics of 
weakly interacting BECs were realized . 

In this letter we shed new light on the interplay be- 
tween disorder and interactions by studying trapped 
BECs under the influence of a disordered potential and a 
one dimensional (ID) optical lattice (OL). The OL cre- 
ates a periodic potential and the randomness of the disor- 
dered potential leads to AL for noninteracting particles 
0. We study how the presence of interactions affects 
nontrivial localization in our necessarily finite system. 

Our experiments were performed with ^^Rb Bose- 
Einstein condensates in an elongated magnetic trap (MT) 



with axial and radial frequencies of = 2t: x 14 Hz and 
iLj± — 271 X 200 Hz, respectively. Further details of our 
experimental apparatus were described previously |l9l |. 
The number of condensed atoms A'^ was varied between 
1.5 • 10^ and 8 • 10^. The OL was provided by a retro- 
refiected laser beam at A = 825 nm superimposed on 
the axial direction of the magnetic trap. The depth of 
the OL was typically set to 6.5 Er, where the recoil en- 
ergy is given by Er = h'^k'^/2m. For this configuration 
the peak chemical potential varied between 0.25 Er and 
0.5 Er- The disorder potential (DP) was produced by 
projecting the image of a randomly structured chrome 
substrate onto the atoms giving rise to a spatially vary- 
ing dipole potential along the axial direction of the cloud. 
Due to the resolution of the imaging system the minimal 
structure size of the DP was limited to 7 /j,m. We define 
the depth of the DP as twice the standard deviation of 
the dipole potential, analogously to 0|. The combined 
potential allowed for the first realization of an ultracold 
disordered lattice gas. 




FIG. 1: Typical absorption images of a BEC with iV = 7 - 10 
released from the combined MT plus DP (left column) and 
MT plus OL plus DP (right column). The second row shows 
the column density and the third row shows the result of a 
ID simulation. The lattice depth was 6.5 Er and the DP had 
a depth of 0.2 Er. 

After the production of the BEC in the MT, we per- 
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FIG. 2: Size of the central peak after 20 ms of ballistic expan- 
sion versus the number of atoms. The clouds were released 
from the following potentials MT (red o), MT plus DP (black 
□ ), MT plus OL (blue <l), MT plus DP plus OL (green t-). 
The lines correspond to a theoretical prediction (see text). 
The lattice depth was 6.5 Er and the DP had a depth of 
0.1 Er. 



formed the following experimental sequence: We first 
ramped up the OL potential over 60 ms, then the DP 
was ramped up over another 60 ms, followed by a hold 
time of 20 ms. Finally all potentials were switched off 
and the atomic density distribution was measured after 
20 ms of ballistic expansion using absorption imaging. 
Alternatively we performed the same experiment with- 
out the OL. 

Figure n shows typical absorption images for the case 
of DP only and for the case of combined DP and OL. 
The obtained density distributions show two character- 
istic features. On one hand they display pronounced 
fringes and on the other hand the axial size of the cen- 
tral peak is modified with respect to the case without 
DP. We extract the axial size of the peak by fitting the 
density with a parabolic distribution. The resulting sizes 
are shown as a function of the atom number in Fig. [3 

Both features can be attributed to the distribution of 
the atoms into the wells of the DP. This can lead to a 
slight fragmentation of the BEG and causes strong fringes 
in the absorption images. These results are in good qual- 
itative agreement with a numerical simulation based on a 
ID Gross-Pitaevskii equation (GPE) as shown in Fig. ^ 
Note that the pronounced interference fringes in the sim- 
ulation result from the interaction dominated axial ex- 
pansion within our ID model. The additional axial con- 
finement due to the DP also leads to an increase of the 
axial size after expansion shown in Fig. [2| The red and 
blue curves show a theoretical prediction based on the 
Thomas-Fermi (TF) approximation. For the black and 
green lines the same functional dependence was fitted to 
the experimental data. This revealed an increase in axial 



size by 25% and 28% respectively. We have used a 3D 
simulation to confirm that this increase is consistent with 
the modification of the chemical potential, introduced by 
the DP. Note, that the change in size depends strongly on 
the exact realization of the disorder. Despite these effects 
of the DP, the computed ground states reveal the absence 
of exponentially localized states (see analysis below) and 
we conclude that the observed localization in the absorp- 
tion images is not caused by quantum interference effects, 
i.e. it does not represent AL. 

In order to understand the experimental results, we 
consider an effective ID model. The BEG spreads over 
more than a hundred wells of the OL, each of the wells 
containing several hundreds of atoms. In this situation 
and for depths of the OL and DP studied here the GPE 
is appropriate |2nl |. In oscillator units corresponding to 
the trap frequency the GPE reads 



g2 ^2 

-y + y + V^,cos2(fcx) 



(1) 

where Vq is the depth of the OL while the DP is rep- 
resented by Vdis(a^). The coupling constant g is chosen 
such that the TF radius equals the axial radius of the 3D 
atomic cloud in the experiment (for the case of = 7T0^ 
presented in Fig. ^ we obtain g — 1500). In all further 
cases we have chosen Vq = 6.5 Er- 

The DP in the experiments changes on a scale much 
larger than the lattice spacing and the condensate healing 
length, / = 1/V8 nna, where n is the condensate density 
and a the atomic scattering length. This suggests the 
applicability of the so-called effective mass analysis |2]| . 
We determine the ground state solution of the stationary 
GPE in the form (l)o{x) = VJ7 f {x)uo{x) , where un{x) is 
the Bloch function corresponding to the ground state of 
the OL potential, f{x) is an envelope function and A/" is a 
constant chosen such that 0o is normalized to unity. This 
leads to an effective GPE where the OL potential is elim- 
inated but the mass of a "particle" and the interaction 
strength become modified. For the experimental parame- 
ters the effective mass is m* = 2.56 and the renormalized 
interaction strength for N = 7 ■ 10"* is g* = 2498. 

Due to the large value of g* we may use the TF ap- 
proximation and obtain the envelope function in the form 
|/(a;)p = [n* — x'^ /2 — Vdisix)]/g* , where ^* is determined 
from the condition J \ f{x)\'^dx = 1. The squared over- 
lap of the obtained 0o with the exact ground state of 
the GPE is 0.999 which implies that the effect of the lat- 
tice potential is reduced to a modification of the coupling 
constant for the TF profile of the combined MT plus DP. 
Thus, similarly to the experiments performed in the ab- 
sence of an OL we observe a fragmentation of the 
BEG induced by the DP but this fragmentation does not 
correspond to Anderson-like localization. 

To enter the localized regime, it is therefore necessary 
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FIG. 3: Superfluid fraction as a function of the coupling con- 
stant g obtained from a ID GPE simulation for a pseudo- 
random potential. Full (open) symbols correspond to a trap 
frequency of 2n x 14 Hz (2tv x 4 Hz). 



to introduce a disorder that changes on a length scale 
comparable to the lattice spacing. Due to the limited 
resolution of the DP imaging optics this poses consid- 
erable experimental difficulties. Alternatively one may 
use a pseudorandom potential obtained with the help of 
two, or more additional optical lattices with incommen- 
surable frequencies |23 |. However, even the realization 
of such a fine scale disorder is not necessarily sufficient 
for the observation of non-trivial localization. Indeed, 
for a solution (/jq of the stationary GPE the nonlinear 
term (7|0o(a;)p may be treated as an additional poten- 
tial. When the atoms accumulate in the wells of the ran- 
dom potential, the nonlinear term in the GPE effectively 
smoothes the potential modulations ^ISj. For typical ex- 
perimental parameters the term g|(/)o(x)p dominates over 
Vdis(a;) and consequently the randomness necessary for 
localization is lost. 

This picture is confirmed by analyzing the dependence 
of the superfluid fraction on the coupling constant g 
shown in Fig. |2| To calculate the superfluid fraction we 
have numerically solved the ID GPE in a box with peri- 
odic boundary conditions in the presence of an OL and a 
pseudorandom potential created by two additional opti- 
cal lattices at 960 nm and 1060 nm with depths of 0.2 E^,.. 
The size of the box was chosen to match the size of the 
atomic cloud in the harmonic potential. The superfluid 
fraction is deflned as fs = 2[Eo{v) ~ Eo{0)]/Nv'^ where 
Eq{v) is the ground state energy when a velocity field v 
is imposed on the system (i.e. we compute the ground 
state solution in the form (I)q{x) cxp{ivx) where (^0(2;) ful- 
fills periodic boundary conditions) |2^ . The superfiuid 
fraction remains large for typical experimental parame- 
ters, indicating the absence of Anderson-like localization. 

To overcome the screening of the disorder potential the 
interaction within the atomic sample has to be reduced. 
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FIG. 4: Ground states of the GPE (note the varying logarith- 
mic scales) for a condensate in the combined potential of the 
MT, OL and pseudorandom potential. The depth of the OL is 
6.5 Er while the depths of the additional lattices forming the 
pseudorandom potential are 0.2 Er- The coupling constants 
g for the panels are: 0.5 (a), 8 (b), 256 (c). Oscillator units 
corresponding to a trap frequency of 27r x 4 Hz are used. 



This can be achieved by reducing the number of atoms, 
lowering the trap frequencies or tuning the scattering 
length via Feshbach resonances. We have performed cal- 
culations for a trap frequency of 27r x 4 Hz and a pseudo- 
random potential equivalent to the one used for Fig. O 
For g = one obtains Anderson-like localization of the 
ground state wavefunction which is characterized by an 
exponential localization |(/)o(a;)p cx exp(— |a;— xqI/Z), with 
the localization length I « 0.027. For such a non- 
interacting system, there might exist several localized 
single particle states with an energy close to the ground 
state. For finite observation times, condensation could 
occur into several of these low energy states, and several 
" small" condensates with different condensate wavefunc- 
tions could coexist. Figure0]suggests that the condensate 
wavefunction becomes a combination of these localized 
states due to nonlinear interactions. 

Increasing g causes the ground state to contain a larger 
number of localization centers. However, the localiza- 
tion length in these cases hardly deviates from the non- 
interacting case. When g is of the order of 500 one can 
no longer distinguish individual localized states and the 
clear signature of non-trivial localization vanishes. This 
is consistent with the appearance of a signiflcant super- 
fluid fraction in Fig. 13 The results shown in Fig. 0J; for 
g = 256 correspond to axial and radial frequencies of 
27r x 4 Hz and 27r x 40 Hz, respectively and N = 10*. 
In this case the simulation shows characteristic features 
of Anderson-like localization while these parameters are 
within experimental reach. The scenario of a cross-over 
from the Anderson to the screening regime, presented 
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FIG. 5: Atomic density after 20 ms of ballistic expansion for 
a condensate prepared initially in the states shown in Fig. |1| 
(left column) and without DP (right column). Oscillator units 
corresponding to a trap frequency of 27r x 4 Hz are used. 



here, is one of the most important results of our analysis. 

Our theoretical investigation also shows that the de- 
tection of the onset of Anderson-like localization using 
a measurement of the density distribution after ballis- 
tic expansion might be difficult. We have calculated the 
atomic density profiles after 20 ms of ballistic expansion 
corresponding to the parameters of Fig. 0] Despite a 
striking difference in the ground state wavefunction, the 
width of the envelope of the zero-momentum peak, which 
is related to the localization length I, does not vary sig- 
nificantly as shown in Fig. [S] In addition Fig. [S]: shows 
that the expansion is dominated by the interaction for ex- 
perimentally accessible values of g within our ID model. 
Future experiments on the detection of localization may 
rather rely on a measurement of the superfluid fraction 
(see |2^) in an accelerated optical lattice. 

In conclusion, we have presented a detailed analysis 
of non-trivial localization for slowly varying potentials 
and in pseudorandom potentials in the presence of in- 
teractions. We have shown the absence of localization 
in the experimental case and explained this effect using 
an effective mass approach. For a truly random poten- 
tial a suppression of Anderson-like localization due to the 
screening by nonlinear interactions was found. An anal- 
ysis for small interactions and a pseudorandom potential 
reveals the characteristic features of Anderson-like local- 
ization. The transition from the localized to the screened 
delocalized regime may be detected via an analysis of the 
superfluid fraction. This work paves the way towards 
the observation of Anderson-like localization in an ex- 
perimentally accessible regime. 
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